Charge Transport and Quantum Phase Transitions in Singlet Superconductor - 
Ferromagnet - Singlet Superconductor Junctions 



o 
o 



(N 



o 
o 



c 

o 

(N 
> 

oo 

(N 
O 

O 

-i— > 



i 

c 

o 
o 



X 



Boris Kastening 1,2 , Dirk K. Morr 3 ' 4 , Lambert Alff 2 , and Karl Bennemann 4 
1 Institut fur Theoretische Physik, Technische Hochschule Aachen, Physikzentrum, 52056 Aachen, Germany 
Institut fur Materialwissenschaft, Technische Universitat Darmstadt, Peters enstrafie 23, 64287 Darmstadt, Germany 
3 Department of Physics, University of Illinois at Chicago, Chicago, Illinois 60607, USA 
4 Institut fur Theoretische Physik, Freie Universitat Berlin, Arnimallee 14, 14195 Berlin, Germany 

(Dated: May 12, 2009) 

We study the Josephson current, Ij, in a junction consisting of two s-wave superconductors 
that are separated by a ferromagnetic barrier possessing a magnetic and non-magnetic scattering 
potential, g and Z, respectively. We discuss the general dependence of Ij on g, Z, and the phase 
difference (f> between the two superconductors. Moreover, we compute the critical current, I c for 
given g and Z, and show that it possesses two lines of non-analyticity in the (g, Z)-plane. We 
identify those regions in the (g, Z)-plane where the Josephson current changes sign with increasing 
temperature without a change in the relative phase between the two superconductors, i.e., without 
a transition between a and n state of the junction. Finally, we show that by changing the relative 
phase <j>, it is possible to tune the junction through a first-order quantum phase transition in which 
the spin polarization of the two superconductors' combined ground state changes from (S z ) — to 
(S z ) = 1/2. 

PACS numbers: 74.50.+r, 74.45.+C, 74.78.-w, 85.25.Cp 



I. INTRODUCTION 



Heterostructures consisting of magnetically active lay- 
ers provide new possibilities for manipulating charge (and 
potentially spin) transport and are hence of great inter- 
est for the field of spin electronics^. Josephson junctions 
consisting of conventional s-wave superconductors and a 
ferromagnetic barrier fall into this category, and their 
study has led to the discovery of a number of fundamen- 
tally new phenomena (for a recent review, see Ref.0 and 
references therein). Among these is the transition from 
a 0-state to a 7r-state in junctions with a metallic ferro- 
magnetic barrier, which is accompanied by a sign change 
(and hence directional change) of the Josephson current. 
This transition signifies an intrinsic phase change of ir be- 
tween the superconductors forming the Josephson junc- 
tion which arises from a temperature dependent decay 
length and oscillation length of the superconducting or- 
der parameter inside the ferromagnetic metal. Such a 
transition was predicted theoretically 3 , and subsequently 
also observed experimentally. It was recently shown, by 
using a phenomenological S-matrix scattering formalism, 
that a sign change of the Josephson current with increas- 
ing temperature can also occur in insulating ferromag- 
netic barriers^. Changing the direction of the Josephson 
current by increasing temperature or by varying the rel- 
ative phase between the superconductors, leads to novel 
types of current switches that possess potential applica- 
tions in quantum information technology 1,7 . Since most 
ferromagnetic barriers do not only possess a magnetic 
scattering potential, but also a non-magnetic one, the 
question naturally arises to what extent the interplay be- 
tween these two types of scattering potentials either al- 
ters the effects discussed above or leads to qualitatively 
new phenomena. 



In this article, we study the Josephson current, Ij, in 
a one-dimensional (ID) Josephson junction consisting of 
two s-wave superconductors and a thin (^-function type) 
ferromagnetic barrier (SFS junction). We start from a 
microscopic Hamiltonian in which the barrier possesses 
a magnetic and non-magnetic scattering potential, de- 
scribed by g and Z, respectively, the former being di- 
rectly proportional to the barrier's magnetization. We 
discuss the general dependence of the charge Josephson 
current on g, Z , and the relative phase, <fi, between the 
two superconductors. In particular, we demonstrate that 
in certain regions of the (<?, Z)-plane, Ij varies contin- 
uously with (f>, while in other regions, and particularly 
around Z = g, Ij exhibits discontinuities. We compute 
the critical current, I c , defined as the maximum Joseph- 
son current for a given g and Z, and show that it pos- 
sesses two lines of non-analytic behavior in the (g, Z)- 
plane. These non-analyticities correspond to discontinu- 
ities in the first and second derivative of I c (with respect 
to g or Z). We show that I c exhibits qualitatively dif- 
ferent dependencies on the scattering strength in differ- 
ent parts of the [g, Z)-plane, which opens the interesting 
(and quite counterintuitive) possibility to increase the 
critical current through the junction by increasing the 
junction's magnetization. Moreover, we identify those re- 
gions of the (g, I?)-plane in which the Josephson current 
changes sign (and thus direction) with increasing temper- 
ature without a change in the relative phase between the 
two superconductors, i.e., without a transition between 
a and n state of the junction. In addition, we find 
that while the total spin Josephson current, I s , flowing 
through the junction is zero, there are two contributions 
to I s , arising from the Andreev and continuum states, re- 
spectively, that are equal in magnitude but possess oppo- 
site signs. We show that if these two contributions can be 
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independently measured, this would open new venues for 
employing the combined spin and charge degrees of free- 
dom in such a junction. Finally, we demonstrate that by 
changing the phase 4> between the superconductors, it is 
possible to tune the junction through a first-order quan- 
tum phase transition in which the spin polarization of the 
superconductors' ground state changes from (S z ) = to 
(S z ) = l/2. 

The theoretical methods used in this study provide di- 
rect insight into the explicit dependence of the Josephson 



current and the Andreev states on the magnetization of 
the junction, and into the interplay between magnetic 
and non-magnetic scattering potentials. 



II. THEORETICAL METHODS 

We take the ID SFS junction to be aligned along the 
z-axis and to be described by the Hamiltonian 



H 



dz 



h 2 d 2 z 

2m 



/i + U(z) 



A(z)v4(zW(z) +H.c. 



M(«)-5^£(*)ffa/J^3(«) 

a,/3 



(!) 



where an d ^v(z) are the fermionic operators that 

create or annihilate a particle with spin a at site z, re- 
spectively. A(z) is the s-wave superconducting gap, and 
U(z) — U~o5(z) describes the non-magnetic (i.e., poten- 
tial) scattering strength of the junction at z = 0. With- 
out loss of generality, we choose the magnetization of 
the junction, M(z), to be parallel to the z-axis, i.e., 
M(z) = Mo(0,0, l)5(z). In order to simplify the no- 
tation, we set g e fiBfJ-o/?i = 1- Moreover, to facilitate the 
discussion of the spin structure of the Josephson current 
and the junction's ground state, we choose the quantiza- 
tion direction to coincide with the direction of the mag- 
netic moment, i.e., with the z-dircction. 

In the following, we use two complementary theoreti- 
cal approaches in order to compute the Josephson current 
through the interface: one follows the Blonder-Tinkham- 
Klapwijk (BTK) approach 8 - (see Sec. Ill A| ). the other one 
starts from the quantum mechanical definition of the cur- 
rent operator and computes its expectation value (see 

sec.nrei). 



A. BTK Ansatz 

At the interface between the two superconductors, two 
Andreev bound states^ with energies E a ^ are formed&i^. 
As was discussed befor o 8 ' 11 ' 12 , and as we will explicitly 
show in Sec. Ill B[ the charge Josephson current flows 
solely through these two bound states, and is hence given 



3=<x,P 



E; 



2k B T 



(2) 



where <f> is the phase difference between the supercon- 
ducting order parameters on the left and right side of 
the junction. In order to compute E a ^, we start from 
Eq. (fTJ) and derive the Bogoliubov-de Gennes (BdG) 



equatio n 12 ! 13 ' 14 ! 15 by introducing the unitary Bogoliubov 
transformation 

^T( z ) =22 u n,a( z ) a n + K,p( z )Pl . ( 3a ) 
n 

*Pl( z ) = ~ U nA Z )Pn + V*,a( Z ) a n ' ( 3b ) 
n 

where the sum runs over all eigenstates of the junction, 
and a n , (3 n are quasi-particle operators in terms of which 
the Hamiltonian, Eq. (JTJ), is diagonal. Defining 



u n,j{z)\ 



(4) 



with j = a,j3, the BdG equation is given by 

HjV„,j(z) = E n ^ nJ (z) . (5) 
Here E n j is the energy of the state $ n j, 



Hi 



-A* -HqtH 



M 



(6) 



where the upper (lower) sign corresponds to j = a(/3) 
and 



h 2 d 2 

H = -^-fi + U 6(z) 
H M =M 5(z) , 

'a z<0, 



A = 



A n e"^ z > 



(7a) 
(7b) 

(7c) 



with Ao chosen real. The case n = corresponds to the 
Andreev bound states, and in what follows, we use Ej as 
a shorthand notation for their energies, E j. 

For the bound states' wave-function on the left and 
right hand sides of the junction, ^o.j t L{z) and ^oj t n(z), 
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respectively, we make the ansatz 



0,j,s 



CO 



s=± 



UQ,j,s,6 
V 0,j,s,5 



„5ikf z 



(8) 



with s = L, R, cl,r — ±1, and kp is the Fermi momen- 
tum. Note that the decay length of the Andreev state, 



Kj 1 = hvp I J Aq — E 2 with vp = hkp/m, itself depends 

on Ej. The solutions of Eq. ([5]) are subject to the bound- 
ary conditions 



*oj,l(0) = *oj,«(0) , 
S.*oj,h(0)-S,*oj,l(0) = 
2m / [7 =fM 



(9a) 



ft 2 



t/ ±M, 



*Oj,ji(0) , (9b) 



where again the upper (lower) sign corresponds to j = 
a(/3). In the limit kp 3> Kj, which holds for supercon- 
ductors with coherence length £ c = vp/A ^> l/k F , the 
solution of the BdG equation yields two Andreev states 
with energies 



E, 



VA + btVA- b 



A ^2 [1 + (.9 + Zf] [\ + {g-Zy\ ' 



(10) 



where the upper (lower) sign corresponds to j = Qt(j3) 
and 

A = (1 + Z' 2 - g 2 ) [cos 2 (0/2) + Z 2 - g 2 ] + 2g 2 , (11a) 
B=y/[l + (g + Z)*] [l + (g-Z) 2 ] 

x [cos 2 (<f>/2) + Z 2 - g 2 ] , (lib) 

with g = mMo/h 2 kp and Z = mUo/h 2 kp. Without loss 
of generality we assume g, Z > from here on. While 
Ep does not change sign as a function of <j> (for Z / 0), 
E a changes sign if < g 2 — Z 2 < 1 (hence for g < Z 
or g > y/l + Z 2 , no sign change of either bound state 
takes place). This sign change, which occurs at a phase 
difference <j>" given by cos 2 (0"/2) = g 2 — Z 2 , indicates 
a first-order phase transition in which the spin polar- 
ization of the superconductors' ground state changes, as 
discussed in more detail in Sec. IIII CI 

For the subsequent discussion, it is necessary to con- 
sider the spin structure of the Andreev states. To this 
end, we compute the local density of states (LDOS), 
N(a,z), (i.e., the local spectral function) for the spin-| 
and spin- 1 components of the Andreev states, which, us- 
ing the Bogoliubov transformation presented in Eq. ([3]), 
are readily obtained as 

N(\,z) = W a {z)\ 2 5{uo~E a ) + \v^{z)\ 2 5{uo+E p ) , 

(12a) 

N(i,z) = \u Az)\ 2 t(^E p ) + \v , a (z)\ 2 S(Lo+E a ) . 

(12b) 

Hence, (for E a rp\ > 0) the wave-function of the a state, 



l^o, a): possesses a particle-like spin-J. component, and a 
hole-like spin-f component, i.e., |$o,a) = \p,l) + |^»t)- 
Similarly for the /3-state, \^o,p) = + \h,i). When 
E a changes sign, the occupation numbers of the spin-f 
and spin-| components of the respective wave-function 
are interchanged. We thus find that the LDOS near the 
junction barrier contains four peaks inside the supercon- 
ducting gap, in agreement with the results obtained in 
Refs. [ft 10. As an important check of our calculations, 
we consider the limit <fr = 0, where the Josephson junc- 
tion is identical to a system, in which a single magnetic 
impurity is located inside a ID s-wave superconductor. 
In this limit, Eq. (fl"0|) yields that only one of the An- 
dreev states exists inside the superconducting gap with 
E a < Aq, while the other Andreev state possesses the 
energy Ep — Aq and is thus part of the continuum (see 
Fig. QJa)). These results are in agreement with those of 
T-matrix and Bogoliubov-de Gennes approaches used in 
the context of impurity scattering in s-wave supercon- 
ductors (for a recent review, see Ref. Il6l). We note that 
in this limit, <j> = 0, our results (and those of Ref. Hit ) for 
E a< p do not agree with the findings in Refs. [see also 
Ref. US. 



B. Quantum Mechanical Current 

Using the operator definition of the quantum mechani- 
cal current, we may resolve the spin-"f and spin-| particle 
currents I-\{z) and I\,(z), which in turn allows to compute 
the charge and spin currents via 



Ij(z) = -e[Ii(z) + h(z)}, 

I s (z) = J[I T (*) -/;(*)]• 



(13a) 
(13b) 



In order to obtain appropriately defined current oper- 
ators if (z) and Ii(z), whose expectation values are the 
currents I-\(z) and Ii(z), we note that the density oper- 
ator of spin-er electrons is 



where 



Pa(z) = Pa(z,z)\ z ,__ 



p a {z,z')=^l{z)^{z'). 



(14) 



(15) 



The quantum mechanical particle current operator cor- 
responding to p a (z) then follows from 



ih 



h{z)= ^-[{d z -d z ,)p„{z,z')] 
2m 



ih 
2m 



{[d z ^t(z)]Mz) ~ ri(z)[dMz)}} , (16) 



and its expectation value 



I <7 (z) = -Im(ipl(z)[dMz)]) 
m 



(17) 
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is the corresponding particle current. In what follows, we 
refer to I a {z) as the "conventional" form of the particle 
current. 

After diagonalizing the Hamiltonian with the Bogoli- 
ubov transformation of Eq. the currents given in 
Eq. (|13[) possess in general contributions from both the 
Andreev bound states and the continuum states. The cal- 
culation of the latter is prohibitively cumbersome when 
using the form of I a given in Eq. (fTT|) . It turns out, 
however, that one can use an alternative formulation to 
evaluate I a by defining a "symmetrized" form of the cur- 
rent operator. Specifically, using the anticommutator 



Hl(z),ip(T'( z ')} = ^aa'S(z - z'), 



(18) 



we can write 



(19) 



We next define a "symmetrized" density via 



p7™(z, z') ee \ [^{z)Mz') + M*WlW] > (20) 

where the second term on the r.h.s. of Eq. (j20|) arises 
from the first term via a particle-hole transformation 
(this fact becomes important when discussing the form of 
the LDOS corresponding to p s J m (z, z), see below). Since 
the ^-function in Eq. (TiH|) does not contribute to the par- 
ticle current operator, it immediately follows that the 
densities defined in Eqs. (fT9"|) and (f2"D|) when inserted into 
Eq. (|16p yield the same spin-er particle current operator. 
We can therefore write I a in a "symmetrized" form as 



Uz) 



ih 

2m 



Kd z -d z ,)?r(z,z')} 



~ {^{z)[d z ^{z)\ + ^{z)[dM(z)\ - h.c.} 



(21) 



One then obtains 



W = —Im{^ a (z)d z ^(z) + M^t(z)) (22) 

as the corresponding "symmetrized" form of the spin-cr 
particle current. While the current thus defined is identi- 
cal to the conventional form given in Eq. (fT7|) , it turns out 
that the calculation of the contributions from continuum 
states using the right hand side of Eq. (f2"2"]) is considerably 
simplified. In what follows, we therefore discuss the con- 
tributions of continuum and Andreev states to the charge 
and spin current using the symmetrized from of the par- 
ticle current, Eq. (f2"2"| . Since the contributions to the 
particle current arising from individual (Andreev or con- 
tinuum) states are different in the conventional and the 
symmetrized form, we distinguish between them by de- 
noting with (without) a tilde the current flowing through 



a specific state within the symmetrized (conventional) 
definition. We then obtain 



7 T(i)( z ) = S-%),n(*). 



(23) 



where 



7 T ,„(z) ee I- a (z) tanh + r n<0 [z) tanh ^L, 



I Un (z) ee II Jz) tanh +I v n Jz) t^h 



with amplitude functions 



(24a) 

@En,a 

(24b) 



T U Z ) = (25) 



h 

2m' 

where f = u,v and j = a, (3, and the sum in Eq. (|23j) 
runs over all states of the system (including Andreev and 
continuum states). 

We begin by computing the contribution to the total 
current arising from the continuum states. Using the 
general ansatz for the form of these states described in 
Appendix El we find that the continuum states with mo- 
mentum ±(fci? ± q) are degenerate and that their energy 
E q is given by (in the limit E q <§; p) 



K = |A | 



h 2 kpq 



(26) 



Diagonalizing the corresponding subspace, one obtains 
three important relations 



^ [I hn (z) + I Un (z)) ocsin2g|z|, 

n(E q ) 

[I 1 ,n(z)-I U n(z)]=0, 

n{E q ) 



(27a) 
(27b) 
(27c) 



n(E q 



n(E q ) 



where the sums run over an orthonormal basis of contin- 
uum states with energy E q . The first relation, Eq. (|27aD , 
implies that given the definition of the charge Josephson 
current, 2j, in Eq. (|13a[) . the contribution to Ij from con- 
tinuum states at the interface is identically zero within 
the symmetrized version of the current. As we show be- 
low, this result also holds when the conventional defini- 
tion of the current operator is used. General arguments 
have been put forward that this result arises since the 
density of continuum states in the presence of a barrier 
is unchanged^. However, away from the interface the 
continuum states carry a non-zero charge current since 
charge conservation requires that the (decaying) current 
through the Andreev states be compensated by a cur- 
rent carried through the continuum states. The second 
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relation, Eq. (|27b) . implies that the contribution of the 
continuum states to the spin Josephson current, as de- 
fined in Eq. (|13bj) . is zero at any position z along the 
junction. Finally, the third relation, Eq. J27c[) . when 
combined with Eq. (|24p . yields that at the barrier, there 
are no contributions from the continuum states to either 
the spin- 1 or spin-j particle current. 

Since the continuum states carry no charge current at 
the barrier (i.e., at z = 0) the total charge Josephson 
current is solely carried by the Andreev states and thus 
given by 



r ./^/./(0) = -e[/ T AS (0) + /f a (0)] 



FAS/ 



(28) 



where / A ? are the currents through the Andreev states 
in the symmetric formulation of the current. These are 
given by the n = term in Eq. (|23j) for which one thus 
has 

I? S (z) = ^ Q (z)tanh^ +/ ^(z)tanh^, (29a) 

Jf (*) = Il p {z) tanh ^ + % a (z) tanh (29b) 

It is straightforward to show that the bound state ansatz 
of Eq. © leads to 



Ilp{z) = Il {z)<xe-^\, 



(30a) 
(30b) 



yielding J AS (z) = lf S {z). Together with Eq. (f27b)l . 
this result implies that within the symmetrized form of 
the particle current, Eq. (|22[) . neither the continuum 
states nor the Andreev bound states carry a spin cur- 
rent. Hence, we obtain that the total spin current Is{z) 
defined in Eq. (|13b[) vanishes at any point along the junc- 
tion, in agreement with the arguments in Ref. 19 . This re- 
sult holds even when the system undergoes a first order 
quantum phase transition in which the spin polarization 
of the junction's ground state changes (see Sec. IIII Cp . 

While the symmetrized form of the Josephson current 
used above allowed for a simpler evaluation of the contri- 
butions arising from the continuum states, any physical 
interpretation of the Josephson current has to be based 
on its conventional definition, given by Eq. (|16p , and the 
form of the spin density in Eq. (p~5|) . In what follows, we 
therefore discuss the form of the charge and spin currents 
within the conventional definition, and compare them 
with the symmetrized results presented above. A con- 
nection between the expressions for J CTi „ and l a>n (which 
represent the currents flowing through state n in the sym- 
metrized and conventional definition, respectively) can 
be made by using the following identities 



J2KM)-Il (z)} = Q, 

n 



(31a) 
(31b) 



where, as in Eq. (|23|) . the sum runs over all states of 
the junction. These identities are derived by applying 
d z — d z i to the anticommutator in Eq. (fT8|) . and subse- 
quently setting a' = a and z' = z, and using the form 
of tpi(z) and tpli 2 ) given in Eq. ([3]). One then finds that 
by subtracting the left hand side of Eq. (|31a|l [Eq. (|31bp ] 
from the right hand side of Eq. (I23p for a =t [a =1], one 
obtains the corresponding expressions for l a ^ n . In partic- 
ular, within the conventional definition, the Josephson 
current through the Andreev bound states is given by 



l^{z) = -Il a {z) 1-tanh 



(3E a 



If(z) 



Given the result for Ij in Eq. 
that 



+ Il (j {z) (l + tanh^) 
(l-tanh^) 
+ ISM fl + tanh^ 



(32a) 



(32b) 



, it immediately follows 



Ij = -e[lf s (p)+lt s (p)] , 



(33) 



implying that also within the conventional definition of 
the particle current, the total charge Josephson current 
at z — is solely carried by the Andreev states with 
no contribution arising from the continuum states. This 
result justifies the use of Eq. @ in the BTK approach of 
Sec. Ill Al for the calculation of the total charge Josephson 
current. Moreover, since the charge Josephson currents 
computed within the BdG approach of Sec. Ill Al should 
be the same as that of the quantum mechanical method 
of Sec. Ill Bl one requires [by combining Eqs. (0), 1(55]). 
(f2"9")l . and (|3"0|) ] that the following important identity be 
satisfied 



dEi 



2^.(0), 



(34) 



with j = a, (3. We have carried out an extensive nu- 
merical survey in the parameter space of our system and 
found this relation to always hold. 

We next consider the form of the spin current in the 
conventional definition of the currents. While the result 
of a vanishing total spin current, that we obtained within 
the symmetrized form, also has to hold within the con- 
ventional framework, we find that the contributions from 
the Andreev and continuum states differ. In particular, 
using the results of Eq. (f3"2"]) . we find that the contribu- 
tion of the Andreev states to the spin current, Ig S , is 
given by 



If{z) = h[ll p {z)-Il a {z)} 



(35) 



which in general does not vanish. However, since the 
total spin current still needs to be zero, it immedi- 
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ately follows that within the conventional definition of 
the particle current, a spin current of equal magnitude, 
but opposite sign to Ig S flows through the continuum 
states. Note, however, that while a spin Joscphson cur- 
rent flows through the continuum states, the charge cur- 
rent through these states is zero, as discussed above. The 
contribution to the spin current provided by the contin- 
uum states thus compensates the spin current through 
the Andreev states and leads to a zero total spin current 
[we return to a discussion of these two (opposite) contri- 
butions to the spin current in Sec. IIII B| . In contrast, in 
the symmetrized version, the spin current through each 
of the continuum states and the Andreev states is exactly 
zero. 



therefore physical and in full agreement with the results 
of Ref. 16 for <f) = 0. It then follows that the physically 
correct interpretation regarding the origin of a zero total 
spin current is that both the Andreev states and the con- 
tinuum states carry a spin current of equal magnitude, 
but opposite sign. We propose that this conclusion can 
be tested by using a spin-resolved STM experiment which 
can distinguish between the LDOS presented in Eq. (JT2J) 
on one hand, and that given in Eqs. (|36[) on the other 
hand. 



III. RESULTS 



In order to understand the different origin of the zero 
spin current in these two frameworks, and to gain insight 
into its physically correct interpretation, we consider the 
local density of states of the Andreev states, N sym (<r, z), 
corresponding to the density p s J m (z) = (p s J m (z,z)), 
which is given by 

iV sym (T, z) = l -\u , a {z)\ 2 - E a ) + 8{u> + E a )] 

+ l\voA z )\ 2 lt(u+Ep) + 8(LU-E l} )] , (36a) 

A sym (|,z) = ^\uoAz)\ 2 [S(Lo-E fj ) + S(Lu+E p )] 

+ ^\voAz)\ 2 [6(u+E a ) +5{u;-E a )] . (36b) 

Here, the second term on the right hand side of each 
equation is obtained from the first one via a particle- 
hole transformation, in agreement with the definition of 
p s J m (z, z') in Eq. (j2TJ)) . As a result, the system now pos- 
sesses two sets of two degenerate bound states, i.e., a 
total of four Andreev bound states, each with a spectral 
weight of 1/2. Within each set, the degenerate bound 
states differ by their spin quantum number: one bound 
state possesses a particle (hole) branch which is spin- 
t (spin- J.), and vice versa for the second bound state. 
Using the relations in Eqs. ([3H| and (|3"^|) . and the def- 
inition of the Josephson current in Eqs. © and (|13ap . 
one immediately finds that such a LDOS leads to the 
expressions for I^f(z) given in Eq. (j2"9"|) . Since the spin 
quantum numbers are opposite between the degenerate 
bound states, it naturally follows that the spin current 
through the Andreev states is zero. Note, however, that 
the LDOS given in Eq. (f3"6"]) is unphysical: it does not 
reflect the symmetry-breaking of the ferromagnetic bar- 
rier since for every spin-"]" branch, there exists a degener- 
ate spin-| branch. Moreover, in the limit <j> — (where 
the barrier represents a single magnetic impurity in a 
ID s-wave superconductor), the symmetrized LDOS of 
Eq. (j3"6"| is in disagreement with that obtained from T- 
matrix and Bogoliubov-de Gennes approaches^. In con- 
trast, the LDOS of Eq. (TT2]) which is based on the con- 
ventional definition of the density operator reflects the 
symmetry-breaking of the ferromagnetic barrier. It is 



A. Charge Transport 

The general dependence of the Andreev state energies 
and the Josephson current on the scattering strength of 
the barrier is presented in Fig. [1] where we plot ±E a ^ 
[Fig. [Ha)] and the resulting Josephson current at T = 
[Fig. [Hb)] as a function of cf> for several values of Z and 
g (the columns in Figs. QJa) and [TJ(b) from left to right 
correspond to g = 0, 4, |, 1, respectively, while the rows 
from bottom to top correspond to Z = 0, ^, |, 1, respec- 
tively.) For a purely non-magnetic barrier [i.e., g = 0, left 
column of Fig. [Ha)], we find in agreement with earlier 
results ^ 12 ' 13 ' 14 ' 20 that the Andreev states are degenerate 
with energies 



(0/2) + Z 2 



cos* 



1 + Z 2 



(37) 



This degeneracy is lifted by a non-zero magnetic scatter- 
ing potential of the junction (i.e., g ^ 0) as shown in the 
three right columns of Fig. [TJa). A qualitatively simi- 
lar result was also found in Refs. l6l lTolll7llT9l . Specifically, 
for a purely magnetic scattering potential of the junction 
[i.e., Z = 0, bottom row of Fig. Ufa)], the energies of the 
Andreev states are given by 



E 



a,f3 



A 



1 + .9 2 



cos (cf>/2) T- 9\jg 2 + sin 2 (0/2) 



(38) 



In agreement with the analytical results presented in 
Eqs. (|3"7) and {55]) we find that with increasing g » Z 
[plots in the lower right corner of Fig. [Tfa)] and Z ^> g 
[plots in the upper left corner of Fig. QJa)] the energies 
of both bound states move toward the gap edge. An 
interesting situation occurs for g = Z [plots along the 
diagonal of Fig. Qla)], since in this case, the effective 
scattering strength for the spin- J. and spm-f electrons is 
V c l s = g + Z = 2Z and Vj s = g - Z = 0, respectively. 
In the unitary scattering limit g — Z 3> 1, we then find 
that the energies of the Andreev states are given by 



E a 



cos 2 (0/2) 
" 2Z 



+ o(z- 3 ) , 



(39a) 



7 




FIG. 1: (Color online), (a) Energies ±E a ^ of the Andreev states. The columns from left to right correspond to g — 0, |, |, 1, 
respectively, while the rows from bottom to top correspond to Z — 0, |, |, 1, respectively. The energies of the a and (3 states 
are indicated by red (grey in print) and black lines, respectively, while the spin-| and spin-j components are indicated by 
dashed-dotted and solid lines. For g = (left column) the a and j3 states are degenerate, (b) The resulting Ij at T = for 
several g and Z as a function of cj>. 
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E{3 

A 



sin 4 (0/2) 
%Z~ 2 



0(Z~ 



(39b) 



suit of the included anomalous Greens functions, one ob- 



tains E a ~ g , as in Eq. 



which immediately leads 



Hence, in the limit g = Z — > oo, the a-state becomes a 
zero energy (midgap) state, while the /3-state moves into 
the continuum. This analytical result is confirmed by the 
numerical results shown in the plots along the diagonal 
of Fig. Ha). 

We next discuss the form of the Josephson current, 
resulting from the form of the Andreev states shown in 
Fig. [2(a). In the unitary scattering limit Z max{j, 1} 
the Josephson current at T = is given by (to leading 
order in Z) 



eAo sin 
h 2Z 2 



(40) 



while for g 3> max{Z, 1} one obtains to leading order in 
9 



Ij 



eAo sin c 



(41) 



Hence, the Josephson currents for a predominantly 
non-magnetic [Eq. (|40|) ] and predominantly magnetic 
[Eq. (|4ip ] barrier differ by a phase shift of 7r in the unitary 
scattering limit. This result also follows from a compar- 
ison of the I j plots in the upper left (Z g) and lower 
right (g 3> Z) corners of Fig.QJb). 

Whether the Josephson current in the junction con- 
sidered here is carried by Cooper pairs, or by single 
electrons, depends on the relative strength of g and Z. 
We first recall that in a purely non-magnetic junction 
(i.e., g = 0) it was argued that the dependence of the 
Josephson current on Z in the unitary scattering limit, 
I j ~ Z~ 2 [see Eq. (|40p], implies that the current is car- 
ried by Cooper pair o 13 > 14 . In contrast, in Josephson junc- 
tions consisting of unconventional superconductors the 
scaling of the Josephson current, Ij ~ Z^ 1 , implies that 
it is carried by single electrons^. Here, we find that for 
a predominantly magnetic junction with g 3> max{Z, 1}, 
Ij also scales with the inverse square of the scattering 
strength (see Eq. (|4"Tj) and the current should thus also 
be carried by Cooper pairs. In contrast, for the case 
Z = g —y oo, we obtain Ij ~ Z~ l and the Joseph- 
son current should thus be carried by single electrons. 
Further support for this conclusion comes from consid- 
ering the dependence of E a on the scattering strength 
in the limit = 0. As mentioned above, in this case, 
the junction is identical to a static impurity in a (ID) 
s-wave superconductor—. If the impurity is purely mag- 
netic, the Andreev state a (which is better known in 
this context as a Shiba state) is formed through scat- 
tering processes involving the creation and destruction 
of Cooper-pairs. This immediately follows from a dia- 
grammatic derivation of the scattering T-matrix which 
includes diagrams that contain the anomalous Greens 
functions, F and F* , representing the creation and de- 
struction of Cooper-pairs, respectively^ 1 -. As a direct re- 



to Ij ~ g for 0^0. In contrast, for an impurity with 
scattering strength Vjg = 2Z and Vj s = 0, the T-matrix 
is given by a series of diagrams that contain the normal 
Greens function only (diagrams containing F and F* are 
forbidden), and hence E a ~ Z^ 1 . Note that for the same 
reason, Josephson junctions consisting of unconventional 
superconductors possess Andreev states whose energies 
also scale as -E 1 , ~ Z^ 1 . This connection between the 
scaling of E a and Ij for 0^0 with the case of impurity 
scattering for = demonstrates that the presence or 
absence of scattering diagrams involving the anomalous 
Greens function determines the nature of the Josephson 
current. 

While the Josephson current is a continuous function 
of for certain combinations of Z and <?, Ij also exhibits 
discontinuities, in particular in the vicinity of Z rs g, 
as shown in Fig. HKb). These discontinuities arise from 
a zero energy crossing of an Andreev state at a certain 
phase, 4>lc, where dE a ^/dcf) ^ 0, as shown in Fig.QJa). 
Since at T = 0, only the negative energy branches of the 
bound states are populated and thus contribute to Ij, 
dE a ^/d(j> leads to a discontinuity in the Josephson 
current at 4>lc- In Fig. [Ha), we present a contour plot 
that represents in which parts of the (g, Z)-plane, Ij ex- 
hibits a continuous or discontinuous dependence on 0. 
In the white regions of Fig.[2](a), Ij is continuous, while 
in the grey and black regions, which are located in the 
vicinity of the Z = g line, it exhibits discontinuities as 
is varied. 

This change between continuous and discontinuous be- 
havior of Ij in the (<?, Z)-plane leads to an interesting 
form of the critical current, I c . Here, we define I c for a 
given g and Z as the maximum absolute value of Ij for 
any < < tt, i.e., I c = max^ [|/j(0)|]. At that value 
of 0, for which |7j(0)| exhibits the maximum value, Ij 
either possesses either a continuous extremum or a dis- 
continuity. In the light (dark) grey areas of Fig.[2(a), I c is 
realized at a continuous extremum (discontinuity) , while 
the discontinuity (continuous extremum) realizes a local 
maximum of \Ij\ only. In the black regions of Fig.fSfa), 
no continuous extremum exists. In order to investigate 
the origin of the qualitatively different behavior of Ij in 
the white, grey and dark regions, we consider six pairs 
of values (g,Z), denoted by the dots in Fig. [2Ja), and 
present the resulting Josephson currents as a function of 
in Fig. [2Kb). We find that the different behavior in the 
grey and black regions arises from a shift in the values of 
at which the discontinuity in Ij occurs. As one moves 
from point (1) to (6) in Fig. [5] two discontinuities first 
emerge at = w, and then move symmetrically toward 
= and = 2-7T, respectively. 

In Fig. (2Ia), we present a contour plot of the critical 
current in the (g, Z)-plane, with white (dark grey) areas 
indicating a large (small) critical current. An analytical 
expression for the critical current in the (g, Z)-plane can 
be obtained along the lines g = and for Z = 0, where 
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FIG. 2: (a) Behavior of Ij((j>) in the (g, Z)-plane. White re- 
gions: I j is continuous. Light (dark) grey regions: Maximum 
of \Ij(4>)\ is realized at a continuous extremum (discontinu- 
ity). Black area: Only extrema at discontinuities exist. The 
white regions are bounded by Z 2 = g 2 and by Z 2 — g 2 — 1. 
(b) Ij(<f>) at T — for the six parameter pairs (g, Z) indicated 
by dots in (a) . 



one finds 




I c (g=0,Z) =— ± 



i c (g,z=o) 



where g m ~ 0.915186 is the largest real solution of the 
equation 4g 6 - Ag 2 + 1 = 0. 

Even though the critical current is a continuous func- 
tion of g and Z, it possesses two lines of non-analyticity 
in the (g, Z)-plane that asymptotically approach Z = g 
for Z, g — > oo. These lines are represented in Fig. [3Ja) as 
solid black lines. Line (1), corresponding to the boundary 
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FIG. 3: (Color online), (a) Contour plot of I c as a function of 
g and Z [white (black) areas indicating a large (small) critical 
current] . The solid black lines indicate non-analyticities in the 
critical current, (b) I c as a function of g for several values of 
Z. (c) I c as a function of Z for several values of g. 



between the light grey and the black region in Fig. [2^a), 
represents a discontinuity in the second derivative of I c . 
In contrast, line (2), which corresponds to the boundary 
between the light grey and dark grey region in Fig. HJa,) 
represents a discontinuity in the first derivative of I c . 
Line (2) also represents a sign change in that value of Ij 
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ducting gap. In order to understand this sign change, 
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FIG. 4: Ic as a function of g for Z — 20. The plot of I c as 
a function of Z for g — 20 is virtually indistinguishable horn 
this plot. 



which determines I c . In other words, to the left (right) 
of line (2), I c is realized by a positive (negative) value 
of I j. These non-analyticities become particularly ap- 
parent when one plots the critical current as a function 
of g (for constant Z) and Z (for constant <?), as shown 
in Figs.[3Jb) and (c), respectively [in Fig. [3fb), we indi- 
cated the non-analyticities for the curve with Z = 1/2 
by arrows]. 

The different dependence of Ij on the scattering 
strength in the limits j > Z or Z > j [see Eqs. (|40]) 
and (|4ip] on one hand and g = Z 1, where 



Ij = 



eAg sin < 



h 4Z ' 
p _ eA sin 0(1 



COS( 



16Z 2 



(44a) 
(44b) 



on the other hand leads to the interesting possibility to 
increase the critical current by increasing the magneti- 
zation, and hence the scattering strength of the bar- 
rier. Specifically, we find I c = ^f[-g~ 2 for g 3> Z and 
I c = ^Z^ 2 for Z > g, while for Z = g > 1, we have 
I c = ^-Z- 1 . As a result, we find that for fixed Z > 1 
and increasing 5, the critical current exhibits a maximum 
at Z = g. This effect is demonstrated in Fig. 0] where we 
present I c as a function of g for Z = 20. Hence, for 
a given non-magnetic scattering strength of a paramag- 
netic barrier, it is possible to increase the critical current 
by increasing the magnetization of the barrier and thus 
its magnetic scattering strength. This increase of the 
magnetization can, for example, be achieved by applying 
a local magnetic field via atomic force microscopy^,. 

We next consider the temperature dependence of the 
Josephson current. Here, we find that the splitting of 
the Andreev states for a non-zero magnetic scattering 
strength can lead to an unconventional temperature de- 
pendence of Ij in which it changes sign with increas- 
ing temperature without a change in the relative phase, 
(j), between the two superconductors. This temperature 
dependence is demonstrated in Fig. G3a), where we as- 
sume a BCS temperature dependence of the supercon- 
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FIG. 5: (Color online), (a) Ij as a function of T/T c for <j> = f, 
Z — 1/2 and several values of g. (b) ±E a: p as a function of 
<f> for Z = 1/2 and g — 0.93. The dotted line corresponds to 
the case <j> — ~ of (a). 

we consider the ^-dependence of E aj p which is shown in 
Fig. [5[b). At T = 0, only the branches indicated by 1 
and 2, belonging to Andreev states f3 and a, respectively, 
are occupied. Since the derivatives dE a jd<\> and dEn/d<f) 
possess opposite signs, the Josephson currents through 
them, IJ < and Ij > 0, flow in opposite directions with 
\Ij \ > \Ij\. Since with increasing temperature, the oc- 
cupation of branches 2 and 3 changes more rapidly than 
those of branches 1 and 4, it follows that the magnitude 
of Ij decreases more quickly than that of IJ. As a result, 



the total current, / 



J 



Ij, eventually changes sign. 
A possible sign change of Ij with increasing tempera- 
ture was previously also discussed in Refs. 5,6. However, 
due to the differences between our results for E a> p and 
those in Refs. [see also Ref. [I?} it is presently unclear, 
whether the origin of the sign change in Refs. 516 is the 
same as the one discussed here. 

The qualitative nature of the temperature dependence 
can be altered via a change of the couplings g and Z, 
as follows from Fig. 02a). We find that in general, a sign 
change of Ij with increasing temperature occurs when (a) 
the particle-like components of both Andreev states pos- 
sess the same spin-polarization and (b) <j) is chosen such 
that the energy difference between the Andreev states 
is sufficiently large. These two conditions can easily be 
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satisfied if at least one of the Andreev states exhibits a 
zero-energy crossing, and <j) is chosen to be close to that 
crossing in the region where (S z ) = 1/2. A zero energy 
crossing, however, occurs only in the grey and black re- 
gions of the {g, Z)-plane shown in Fig. Ufa). Thus, in 
order to observe a sign change of Ij with temperature, 
one should select a barrier whose scattering potentials 
are close to the Z = g line. Note that a similar tempera- 
ture dependent sign change is also predicted to occur in 
Josephson junctions consisting of triplet superconductors 
and a ferromagnetic barrier—. It is important to stress 
that the temperature dependent sign change discussed 
above is qualitatively different from the one reported by 
Ryazanov et al. A . There, the sign change arises from a 
transition of the junction from a 0-phase state at high 
temperatures to a 7r-phase state at low temperatures due 
to a temperature dependent coherence length 2 . In con- 
trast, in our case the sign change arises from a change 
in the population of the Andreev states, with the rela- 
tive phase between the superconductors remaining un- 
changed. 



B. Spin Transport through Andreev States 

We argued in Sec. Ill Bl that the origin of the zero total 
spin current lies in the fact that the spin current through 
the Andreev states is compensated by a spin current 
through the continuum states that is equal in magni- 
tude, but opposite in sign. The question thus naturally 
arises if it is possible to separately measure these two 
contributions to the total spin current. While one could 
envision several experimental set-ups in which this could 
be achieved, for example, by using two Josephson tunnel- 
ing STM tips, one on each side of the junction, we cannot 
provide a definite answer to this question at the moment. 
However, the ability to measure these contributions sep- 
arately would open new venues for using the combined 
spin and charge degrees of freedom in such a junction. 
In particular, it would be possible to make use either of 
a spin polarized (non-zero) charge Josephson current by 
considering the current through the Andreev states, or 
of a spin Josephson current without a charge Josephson 
current by considering the current through the contin- 
uum states. To exemplify these possibilities, we consider 
in what follows the spin polarization of the Josephson 
current through the Andreev states. We first define a 
spin polarization V of the Josephson current via 




For V = — 1(+1), the Josephson current through the An- 
dreev states is completely spin polarized and thus solely 
carried by spin- J, (spin-|) electrons. In Fig. [6] we present 
V as a function of cf> for g = 1/3, Z = 2/3 and T = (the 
corresponding charge Josephson current is shown in the 
second panel of the second row in Fig. [TJ)). We find that 
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FIG. 6: Spin polarization V of the Josephson current through 
the Andreev states at T — as a function of <f> for g — 1/3 
and Z = 2/3. 



V is non-zero for all 4>, and its magnitude reaches a maxi- 
mum at <f> = 7r /2 (note that while V — > — 1 at — > 0, one 
has at the same time Ij ->• 0). At T = 0, lf s {lf s ) is 
carried solely by the occupied particle-like branch of the 
/3-state (a-state), as follows immediately from Eq. (|32[) . 
This is consistent with the observation in Eq. (fT2]) that 
the particle-like component of the a and /3 states pos- 
sess the spin quantum number S z = —1/2 (spin- J, ) and 
S z = +1/2 (spin-t), respectively. Since \I^ a \ > |7^| for 
all <fr, one finds that the Josephson current through the 
Andreev states is partially spin-| polarized, as shown in 
Fig. [S] The degree of spin-polarization varies with <f> due 
to the changing relative contributions of Iq a and I£ p 
to Ij. Finally, for g = 0, the two Andreev states are 
degenerate, and hence V = 0. 



C. First Order Quantum Phase Transition 

Another interesting effect arising from the combina- 
tion of magnetic and non-magnetic scattering strength 
of the barrier is the possibility to tune the Josephson 
junction through a first-order quantum phase transition 
in which the ground state of the entire system (i.e., the 
combined ground state of both superconductors) changes 
its spin polarization from (S z ) = to (S z ) = 1/2 (assum- 
ing without loss of generality that the magnetization of 
the barrier points along the £-axis) . This type of first or- 
der transition, which is well known from static magnetic 
impurities in s-wave superconductors, where it was first 
discussed by Sakurai2i, can occur either with increasing 
scattering strength, </, of a magnetic impurity 2 ^, or due 
to quantum interference effects 2 ^; its fingerprint is a zero- 
energy crossing of the impurity induced fermionic states 
inside the superconducting gap. The phase transition 
arises from a level crossing in the superconductor's free 
energy, J 7 , resulting in a discontinuity of d? ' /dJ at the 
transition; hence the first order nature of the transition 
(for a more detailed discussion see Ref.OJ). The Joseph- 
son junction considered here provides a new possibility 
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to tune the system through such a first-order phase tran- 
sition by varying the phase difference between the super- 
conductors. In Appendix [HI we explicitly show that the 
change in the spin polarization of the junction coincides 
with the zero energy crossing of the a state, and is thus 
solely driven by the Andreev states. In contrast, the con- 
tribution of the continuum states to the spin polarization 
vanishes. 

As an example of such a phase transition we consider 
the case Z = 1/3 and g = 2/3, for which the energies of 
the Andreev states and the resulting Josephson current 
are shown in the third panel in the third row of Figs. [Ha) 
and (b), respectively [we indicate in the panel in Fig.QJa) 
the spin quantum number of all the components of the 
Andreev states]. With increasing <f>, the a state crosses 
zero energy at <p% — 2 arccos \J g 2 — Z 2 0.6087T, such 
that for <fi > <j>", the spin-| component of the a state 
is particle-like, while its spin-| component is hole-like. 
This transition results in a change of the ground state 
spin polarization from (S z ) = to (S z ) — 1/2. More- 
over, the a state crosses zero-energy again at <f>"' — 
2[tt - arccos V ' g 2 - Z 2 ] « 1.392tt, such that for 4> > 
its spin-"f component is hole-like, while its spin- J, compo- 
nent is particle-like. As a result, the spin ground state 
changes from (S z ) = 1/2 back to (S z ) = 0. The range 
of <fr for which (S z ) = 1/2 is indicated in the panel of 
Fig. [Hb) by dotted lines. In general, one finds that 
(S z ) = 1/2 for those (j) which satisfy cos 2 (</>/2) < g 2 -Z 2 . 
Consequently, as g is further increased (keeping Z fixed) , 
the range of <f> for which (S z ) = 1/2 increases. When 
g exceeds the upper critical value g^ = \/l + Z 2 , one 
finds (S z ) = 1/2 for all <f>. In contrast, when g is smaller 
than the lower critical value gf = Z, one has (S z ) = 
for all (j). As already mentioned above, the total spin 
current through the junction is zero even in a state with 
(S z ) = 1/2. 

IV. CONCLUSIONS 

In summary, we have studied the Josephson current, 
Ij, in a ID Josephson junction consisting of two s-wave 
superconductors and a thin ((S-function type) ferromag- 
netic barrier. To this end, we used two complemen- 
tary theoretical approaches: the BTK method, and an 
approach starting from the quantum mechanical defini- 
tion of the current operator. We discussed the general 
dependence of the charge Josephson current on g, Z, 
and the relative phase, <j), between the two superconduc- 
tors. Specifically, we showed that in certain regions of 
the (g, Z)-plane, Ij varies continuously with </>, while in 
other regions, and particularly around Z — g, Ij exhibits 
discontinuities. We computed the critical current, I c , 
defined as the maximum Josephson current for a given 
g and Z, and we showed that it possesses two lines of 
non-analytic behavior in the (g, Z)-plane. These non- 
analyticities correspond to discontinuities in the first and 
second derivative of I c (with respect to g or Z). We 



demonstrated that I c exhibits qualitatively different de- 
pendencies on the scattering strength in different parts 
of the (g, Z)-plane, which opens the interesting possibil- 
ities to increase the critical current through the junction 
by increasing the junction's magnetization. This effect 
possesses potential applications in the fields of quantum 
information technology^. Moreover, we showed that for 
certain values of g, Z, the Josephson current changes 
sign (and thus direction) with increasing temperature 
without a change in the relative phase between the two 
superconductors, i.e., without a transition between a 
and 7r state of the junction. We showed that this sign 
change is entirely due to a temperature-dependent change 
in the occupation of the Andreev states. In agreement 
with earlier results, we demonstrated that the contin- 
uum states do not contribute to the charge current, and 
that therefore the charge Josephson current is carried en- 
tirely by the Andreev states. We also showed that while 
the total spin Josephson current through the junction 
is zero, the Andreev states and the continuum states 
separately carry a non-zero spin current of equal mag- 
nitude but opposite sign. The possibility of measuring 
these contributions separately would open new venues 
for employing the combined spin and charge degrees of 
freedom in such a junction for potential applications in 
spin electronics and quantum information technology. Fi- 
nally, we demonstrate that by changing the phase <j) be- 
tween the superconductors, it is possible to tune the 
junction through a first-order quantum phase transition 
in which the spin polarization of the superconductors' 
ground state changes between (S z ) = and (S z ) = 1/2. 

Experimentally, the effects discussed here could be 
studied in junctions with magnetically doped insulating 
barriers based on MgO, ZnO, or TiC>2- In these materials 
one can imagine to vary g independently via the substitu- 
tion of magnetic dopants such as Co, Mn, etc. and/or by 
changing their concentration, or by applying a small mag- 
netic field, for example, via atomic force microscopy^. 
Moreover, Z can be altered by the choice of material 
and the junction width. It is possible to control the bar- 
rier width of complex oxides using layer-by-layer growth 
techniques monitored by reflection high energy electron 
diffraction (RHEED) on the unit-cell level, which is much 
smaller than the coherence length of a typical s-wave su- 
perconductor. Hence, we expect that the results derived 
above for a (5-functional barrier should be observable in 
experimental systems with a non-zero barrier width, d, 
as long as d is much smaller than the superconducting 
coherence length. Note that the novel Josephson behav- 
ior described in this paper will occur in addition to the 
effects that are expected from the proximity induced sign 
change in the superconducting order parameter as a func- 
tion of the ferromagnetic barrier thickness^. 

Finally, scattering off the barrier leads to a suppression 
of the superconducting order parameter near the barrier, 
which was not accounted for in the approach presented 
above. However, in s-wave superconductors, the length 
scale over which the order parameter recovers its bulk 
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value near scattering centers (such as a junction) is set 
by 1/kp 27 . This length scale is in general much shorter 
than both the superconducting coherence length £ c and 
the decay length of the Andreev bound states, « , with 



£c < k 



and we thus expect that the inclusion 



of a spatially varying order parameter does not alter the 
qualitative nature of our results presented above. 
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APPENDIX A: CONTRIBUTION OF THE 
CONTINUUM STATES TO Ij AND Is 

For the scattering states' wave-function of Eq. (0| we 
make the ansatz 



E 

<5,e=± 



u n,j,s,8,e j p i(SkF+tq)z 



(Al) 



with s = L,R referring to the left (z < 0) and right 
{z > 0) hand side of the junction, respectively, with 
j — a, P, with q > 0, and where kp is the Fermi momen- 
tum. The corresponding solutions of Eq. ([5]) are subject 
to the boundary conditions ©. For a given E q > |A|, 
there are eight continuum states with a positive- and a 



negative-energy branch for each. For energies small com- 
pared to the Fermi energy, q is given by (|26p . 



Define 



cosh 7 = — — . 
Ao 



(A2) 



To be specific, consider the negative-energy branches 
with energy — E q . The BdG equations imply 

y w , e = e s ^/\, sA e = eV'-Wvj,.^ , (A3) 

with j = a, (3 and s = L, R and where we have omitted 
the index n labeling the different states of energy —E q . 



With tpi, = and 4>r — —<j) and defining 
e ,r ± _ l+i{Z±g) 

so that 



and 



7T 7T 
77 < T ± < 7T 

2 2 



Z = — (tan r + + tan r_ ) 
1 , 

g = — (tan t + — tan r_ J 



(A4) 
(A5) 

(A6a) 
(A6b) 



we may organize the boundary conditions ^ in matrix 
form by writing 



(A7) 



with 



vis = (yj,s,+,+,yj,s,+,-,yj,s,-,+,yj,8-,-) (A8) 



and 



M, 



V 



e -7/2 

e -i0/2+7/2 
e -ir T -7/2 
-ir± -10/2+7/2 



e 

e" 
e 



+7/2 
i0/2-7/2 
-ir T +7/2 
i0/2- 7 /2 



r 



,+7/2 



g-t^/2-7/2 
_ e ir T + 7 /2 

■e" 



Jr± -10/2-7/2 



e -7/2 
e -i0/2+ 7 /2 
_gir T -7/2 
Jt± -10/2+7/2 



(A9) 



where the upper (lower) sign corresponds to a (/?). Thus 
the negative-energy scattering states for a given energy 
— E q may be represented by yj. Reinstating labels m 
and n for the different states for a given energy, we de- 
mand that for continuum-normalized orthogonal states 



y m< j and y n j holds 

ym,j,s,S,eyn,j,s,5,e = z \n^Q j z n.j 

s=LM8,t=± 

ftrnnCr^Eq} , 



(A10) 
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with 



Qj = (MM)- 1 + (MjMh- 1 * 



(All) 



and where C r (E q ) is real and may be chosen to depend 
on the energy E q . That is, finding an orthonormal ba- 
sis of scattering states for a given energy boils down to 
diagonalizing the Qj and using (|A3|) and (1A7|) to obtain 
all corresponding coefficients u n _j tS ^^ and v n ,j, s ,s,e- Nu- 
merically, we find that 



(\ U n,j,s,+,±\ 

n(E q ) 



L n,j,s, — ,=j= | 



l ) = 0, (A12a) 



and the same for the v n j tS § <e . Moreover, we numerically 
obtain 



*n,j,s, — ,-f- 



n(E q 



= djA r (E q ) + ic s Ai(E q ), 

(A13) 

and the same for the V nt j <s ,s,e, with d a ^p = ±1, cl.r = 
±1, and where A r (E q ) and Ai(E q ) are dimensionless and 
generally nonzero real coefficient functions. 
It is straightforward then to show that 



Im/u* a d z u n , a ± v* n a d z v n ^ a 

n(E q ) 

± u* n [j d z u n ^ + v* tj3 d z Vn t p 

__ [%kAi(E q )sva2q\z\, 
~ 10, 



(A14) 



leading immediately to the results (J27 



APPENDIX B: SPIN GROUND STATE OF THE 
JUNCTION AND FIRST ORDER PHASE 
TRANSITION 



with p a (z,z') from (JT5J) . It is then straightforward to 
show that 



Ps(z) = (ps(z)) 

_ h 

~ 2 



^ S [Pn,«,aW + P7i,t;,a(z)] tanh 



[Pn,t>,/9(-?) + Pn,v,,p(z)] tanh — ^- (B2) 



where 



(B3) 



0, (A12b) witn f = u ,v and j = a, (3. 



Note that for any given energy, there exist pairs of con- 
tinuum states that differ only in the spin quantum num- 
bers of their hole-like and particle-like branches. This re- 
sult, combined with the normalization of the continuum 
states then yields that their contribution to the spin of 
the junction vanishes. Using next the normalization of 
the bound states 



+ OC 



dz[ulAz)u ,j{z) + vlJz)v 0d (z)} = 1, (B4) 



with j = a, (3, we obtain that the spin of the system is 
solely determined by the Andreev bound states and given 
by 



+ OC 



dzpf{z) 



- — (tanh — — — — tanh — — — | 
4 V 2 2 J 



(B5) 



Since our conventions are such that always Ep > 0, but 
the sign of E a can vary, we obtain for T = 



The total spin of the system receives contributions 
from both the Andreev bound states and the scattering 
states. Define the spin density by 



Ps(z)= -[p](z,z') -pi(z, z')} 



(Bl) 







E a > 0, 



h/2 E a < 0, 



(B6) 



signaling a quantum phase transition caused by the zero- 
energy crossing of one of the Andreev bound states. 
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